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We evaluate the charge and spin susceptibilities of the 2D attractive Hubbard model and we 
compare our results with Montecarlo simulations on the same model. We discuss the possibility 
to include topological Kosterlitz-Thouless superconducting fluctuations in a standard perturbative 
approach substituting in the fluctuation propagator the Ginzburg-Landau correlation length with 
the Kosterlitz-Thouless correlation length. 
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The discovery of spin and charge pseudogaps in the 
normal state of underdoped superconducting cuprates |j]] 
has triggered a renewed interest on the physics of pre- 
formed Cooper pairs. The actual source of the pseudo- 
gaps (pairing, and/or spin-, and/or charge fluctuations) 
and the leading mechanisms responsible for the reduc- 
tion of the superfluid density at low temperature (clas- 
sical phase fluctuations, collective modes, quasiparticle 
excitations) are still debated. However, many indica- 
tions support the idea that pairing occurs below some 
crossover temperature T*, while the phase coherence is 
established at a sizable lower temperature. The low den- 
sity of carriers resulting in a low superfluid density and 
the short coherence length £ ~ 10-^20A, support the rel- 
evance of the superconducting phase fluctuations in the 
thermodynamic and dynamic properties of these mate- 
rials. Moreover, although no discontinuity of the super- 
fluid density at T c is observed, the strong anisotropy of 
the cuprates suggests that some features of a Kosterlitz- 
Thouless (KT) transition could be present in these sys- 
tems g]. Therefore it is worth investigating the effects 
of the topological vortex-antivortex phase fluctuations on 
the various properties of a 2D superconductor. In partic- 
ular, an important issue concerns the inclusion of these 
effects in evaluating thermodynamic quantities like the 
spin susceptibility or the charge susceptibility. In this 
context, the aim of the present work is to look for possible 
connections between the perturbative scheme leading to 
the standard time-dependent Ginzburg-Landau (TDGL) 
results and the KT physics. 

Halperin and Nelson S have shown that, in the KT 
regime, the contributions of superconducting fluctua- 
tions to the conductivity above Tkt have the same func- 
tional form, in terms of the correlation length £, as the 
Aslamazov-Larkin contributions of the standard TDGL 
theory, (Jkt{0 — &gl(£.) ~ £ 2 - The same holds for the 
fluctuation contribution to the diamagnetism Xkt(0 — 
Xgl(0 ~ £ 2 - I n s phe of the same correlation length 
dependence, conductivity and diamagnetism in KT or 



TDGL theory have a completely different temperature 
dependence, induced by the different temperature de- 
pendence of the correlation length in the two theories. 
The KT correlation length diverges exponentially at Tkt 
while the GL correlation length diverges as a power-law 
with the classical exponent v = |. Therefore the KT 
conductivity and the diamagnetic susceptibility diverge 
exponentially at Tkt while the same quantities in the 
TDGL theory diverge as a power-law at T c with a criti- 
cal exponent 7 = 1. In the present work we investigate 
the possibility that, in analogy with conductivity and dia- 
magnetism, the correct behavior of the spin and charge 
susceptibilities in the KT regime can be simply recovered 
by inserting the KT correlation length in their TDGL ex- 
pressions. We shall find that this prescription does work 
for the spin susceptibility while it does not for the charge 
susceptibility. 

We analyze the two-dimensional negative-C/ Hubbard 
model [In] which is the simplest minimal model where 
the distinct occurrence of pairing and phase coherence 
can be investigated. Within this model, the spin suscep- 
tibility Xs and the charge compressibility \c ar e calcu- 
lated on a two-dimensional square lattice by performing a 
loop expansion with the fermions exchanging the Cooper- 
fluctuations propagator in the standard form. Before giv- 
ing the technical details of our treatment, we immediately 
present our results. 

Figure 1 shows the behavior of the spin susceptibility 
when the correlation length is assumed either of the GL 
form (dashed line with crosses) or of the KT form (dot- 
ted line with stars). Both curves are compared with the 
Montecarlo data obtained in Ref. 0W for the negative-17 
Hubbard model with U = — At (t is the nearest-neighbor 
hopping) at filling n — 0.5 electrons per cell. The critical 
temperature Tkt of the KT superconducting transition, 
as extracted from numerical calculations, is Tkt = 0.054 
and has been used as the input critical temperature for 
our perturbative calculations. In the Montecarlo data, 
for T less then T* ~ t ^> Tkt, Xs starts decreasing. This 



indicates the existence of strong superconducting fluctu- 
ations in the temperature range between the mean-field 
transition temperature (Tbcs — 0.6i) and the true KT 
transition. It is apparent from Fig. 1 that the rapid de- 
crease of the spin susceptibility in the Montecarlo results 
is well fitted by inserting in the correlation length the KT 
temperature dependence as given by the expression 



£,kt{T) =£ c exp 
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Here £ c is an effective size of the core of the vortex that 
we take of the order of the zero temperature correlation 
length £oi and b is a positive constant of the order of 
unity. This specific form of the KT correlation length 
has been derived along the line of Ref. 0, although it 
differs slightly from the one commonly quoted in the lit- 
erature H||. We shall comment on this later. Notice 
that the KT mass term (inverse square of the correla- 
tion length) of the Cooper propagator remains small and 
generates strong fluctuations, in a wider range of temper- 
atures than the GL mass with the same critical tempera- 
ture in agreement with Montecarlo data. The GL correla- 
tion length is instead completely inadequate to reproduce 
the Montecarlo data in the all range of temperatures. 



FIG. 1. Comparison between Montecarlo spin susceptibil- 
ity (taken from Ref. [6]) and the spin susceptibility cal- 
culated using the Ginzburg-Landau (X? ) and the Koster- 
litz-Thouless (x s T ) correlation length. 

The fit in Fig. 1 stops at T ~ O.lt because there are 
no numerical data below this value. This also appears to 
be the lower limit for our approach to work. Indeed for 
T ~ 0.09t the TDGL expression for Xs develops a non 
physical behavior (x s < 0), indicating that the perturba- 
tive scheme no longer applies near Tkt- Whit this cau- 
tion in mind, the results of Fig. I indicate that the simple 
loop expansion we adopted is able to reproduce the spin 
susceptibility in a wide range of temperatures. They sup- 
port the idea that the main effect of the vortex-antivortex 



phase fluctuations on the spin susceptibility is embedded 
in (and satisfactorily accounted for by) the temperature 
dependence of the £kt{T) correlation length, in analogy 
with the conductivity and diamagnetism. 

On the other hand, as seen in Fig. 2, the same type of 
calculations for the charge susceptibility fail in describ- 
ing the nearly constant (but with sizeable error bars) be- 
havior obtained numerically. In particular, we find that 
the Aslamazov-Larkin (AL) contribution, which does not 
contribute to the spin susceptibility, strongly enhances 
X C (T) and eventually leads to a divergent \c near Tkt- 
As a consequence Xc(T) strongly deviates from the Mon- 
tecarlo results for T < Tbcs- In Fig. 2 we also report the 
RPA resummation of the bare bubble in the charge chan- 
nel that fits the available Montecarlo data, to obtain, by 
extrapolation, the Xc(T) at higher temperature. 
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FIG. 2. Comparison between Montecarlo charge suscep- 
tibility (taken from Ref. [8]), the charge susceptibility cal- 
culated using the Kosterlitz-Thouless correlation length (x c ) 
and the RPA resummation of the bare bubble (x? P )■ 



With this respect Xc appears to behave as the specific 
heat c v , for which the 2D-TDGL expression c v ~ £, GL 
does not reproduce the correct KT result c v ~ £kt> 
even when expressed in terms of the correlation length. 
For the specific heat this happens despite the free en- 
ergies in the two theories have the same leading behav- 
ior when written in terms of the respective correlation 
lengths: F GL ~ £, G l Infer, and F KT - Q 2 T §,|. Indeed, 
since c v involves the second derivative of F with respect 
to temperature the different temperature dependences of 
the correlation lengths (and the subleading In ^ql factor) 
lead to completely different results in the two theories. 
Our result for Xc has the same origin: The charge re- 
sponse at w = 0, q — > can be obtained as a chemical 
potential derivative of the free energy. Now, since the 
critical temperature depends on the chemical potential 
T c = T c {(x), a total derivative with respect to \i also in- 
volves derivatives with respect to T Cl and, in turn, deriva- 
tives of £. Therefore the temperature dependence in Xc 



not only arises from the temperature dependence of £(T), 
but also depends on d^/dT c . In fact one gets the same 
TDGL singular contribution ~ £ 2 for Xc and c v . Our 
simple perturbative expansion, where the leading tem- 
perature dependence only arises from the mass term £ -2 
of the Cooper fluctuation propagator in the TDGL ex- 
pression, fails to reproduce the correct temperature de- 
pendence for \c in the same way as it fails in evaluating 
the specific heat. 

We now describe the details of our calculations. The 
model we consider is given by 
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where t is the hopping between nearest-neighbor sites, 
U < the strength of the attraction and /i the 
chemical potential. The standard ladder resumma- 
tion of diagrams leads to the Cooper pair propagator 
L(q, il t ) = -U/ (1 + Uxl P {q, Hi)) where y^(q, Hi) is the 
bare particle-particle bubble, being q the momenta and 
Qi the Matsubara frequency. In the normal state, within 
the standard GL approach, at small q and Qi one has 



L- 1 (q,n l )^N Q (e + r ] q 2 + 1 \n l |) 



(3) 



where N is the density of states at the Fermi energy, 
77 = 7C(3)/(32^ 2 )(« F /T C ) 2 ~ £ 2 , and 7 = tt/(8T c ). The 
mass term e = ln(T/T c ) = (£/£o)~ 2 of the propagator 
controls the distance from the superconducting transi- 
tion. In the standard GL approach e ~ £q L and near T c 
it goes to zero as (T — T c )/T c . 

We study the charge and spin susceptibilities by eval- 
uating the one loop corrections A^ c (charge channel) 
and Axs (spin channel) to the bare particle-hole bub- 
ble Xo h , X p c h s = X ph + &Xc, s - The charge (c) and spin 
(s) bubbles x ph s are then inserted in the RPA resum- 
mation to get the charge and spin susceptibilities (see 
below). In the one loop expansion, we include diagrams 
containing only one integration on the bosonic variables 
(q, ili) (i.e. one bosonic loop) of the fluctuation propa- 
gator L(q,£li), obtaining three kinds of diagrams which 
contribute differently to the spin and charge susceptibil- 
ities: the selfenergy diagrams, where L(q,£li) renormal- 
izes the one particle bare Green function (DOS contribu- 
tion); the vertex diagrams, where L(q,fli) renormalizes 
the vertex, connecting two bare Green function (Maki- 
Thompson (MT) contribution); the Aslamazov-Larkin 
(AL) diagrams, containing two fluctuation propagators. 
Moreover it is necessary to add the counterterms (CT) 
proportional to the shift of the chemical potential 5/i, 
which is required to preserve the number of particles. We 
notice that the one loop expansion for the charge and the 
spin susceptibilities satisfies the relation, derived from 
spin and charge conservation, Xs,c(q = 0, il =/= 0) = 0. 
One obtains: 



A Xs = 4 • DOS - 2 • MT + 4-CT (4) 

A Xc = 4 • DOS + 2 • MT + 4 • AL + A-CT. (5) 

The absence of the AL contribution and the (oppo- 
site) sign of the MT diagrams in the spin susceptibility 
is the consequence of the vertex spin structure, as shown 
in Ref. DTl. Moreover the leading DOS contributions to 
the charge susceptibility cancel the MT ones. The AL 
diagrams give therefore the most important contribution 
to the charge susceptibility (being the CT diagrams sub- 
dominant respect to them) fljjfl . 

According to the physical assumption outlined above 
that the TDGL and KT temperature dependencies are 
essentially ruled by the correlation lengths, we have al- 
ternatively taken Eq.(g) with £ = £gl and £ = £,kt- 
In the calculation with £gl we used T c = Tkt and the 
mass term e = ln(T/T c ), while in the calculation with 
£,KT we used Eq.(|l|) with b = 1.6 and £ c = £ - I n both 
cases we took the coefficients ij and 7 given by the cor- 
responding expressions reported below Eq.d3) calculated 
with T c = Tbcs- This choice was motivated by the plau- 
sible assumption that 77 and 7 change little once the fluc- 
tuations are predominantly in the phase sector. In any 
case we checked that our results are rather stable with 
respect to modifications of r) and 7. 

The charge and the spin susceptibilities are finally ob- 
tained by the RPA resummation of the corrected charge 

and spin bubbles Xc,s — X ph s l \^ ± (U C}S /2)x P l h s ) where 
the plus (minus) sign is associated to the charge (spin) 
susceptibility. Notice that, following the analysis of Ref. 
M , the RPA expressions of both susceptibilities contain 
an effective local interaction C/ C , s instead of the bare U 
in order to properly fit the high temperature region of 
the Montecarlo data. The validity of the RPA form for 
the spin susceptibility is also found in the context of the 
positive-t/ Hubbard model pa. However, while in Ref. 
H the bare bubbles were resummed and a value U = 6.5 
was obtained for U — —At and < n >— 0.5, in our case 
we resum the bubbles already containing the Ax s correc- 
tions and a different value U s = —4.6 is needed to match 
the RPA calculation with the high temperature Monte- 
carlo data. For the charge susceptibility the comparison 
with the RPA resummation in terms of the x P reported 
in Fig. 2 gives U c = —1.6. 

We now comment on the expression in Eq.l that we 
used for the KT correlation lenght. We wrote this 
expression following Halperin and Nelson J3|. They 
introduce into the KT correlation length £kt — 
aexp \b(nJ /IzbT — 1)] for the classical XY model (with 
coupling J and lattice spacing a) a temperature depen- 
dent J(T) = n s (T)/8m and take a = £ c . Here the super- 
fluid density n s (T) is taken to vanish linearly at a tem- 
perature Tq(> Tkt) to be determined selfconsistently by 
the request that Tq should include the effect of the fluc- 
tuations at scale lower than £ c . Our expression (1) is 



obtained by taking To ~ Tbcs and £,c — Co with the idea 
that phase fluctuations are the most important effect all 
over the range of temperatures Tkt ^ T < Tbcs (at 
least in evaluating Xs and Xc) M l- 

The results of the simple procedure outlined above are 
quite satisfactory for the spin susceptibility. This in- 
dicates that the main temperature dependence of this 
quantity actually arises from the specific KT temperature 
dependence of the correlation length, which thus brings 
along the physics of the vortex-antivortex phase fluctua- 
tions into a simple perturbative scheme. The same is not 
true for the compressibility, as for the specific heat, since 
these quantities also involve temperature derivatives of 
£kt- 

Our method, developed for the 2D attractive Hubbard 
model, can be useful to understand the role of the super- 
conducting phase fluctuations in quasi-2D cuprate super- 
conductors. In this context the recent finding that KT 
signatures, which are absent in the static conductivity, 
are progressively more evident in the dynamical conduc- 
tivity at shorter timescales [|| encourages to extend our 
analisys to other frequency-dependent quantities. In par- 
ticular it is of obvious interest to explore the possibility 
to include in a simple perturbative scheme along the lines 
followed in the present work the effects of KT topolog- 
ical phase fluctuations on dynamical quantities like the 
optical conductivity and single-particle spectra. 
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